Fluid accumulation in thin-film fiows driven by surface 
tension and gravity (I): Rigorous analysis of a drainage 

equation. 
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Abstract 

We derive a boundary layer equation describing accumulation regions within a thin-film 
approximation framework where gravity and surface tension balance. As part of the analysis 
of this problem we investigate in detail and rigorously the following drainage equation 
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In particular, we prove that all solutions that do not satisfy $ ^- 1 as r — >■ oo are oscillatory, 
and that they oscillate in a very specific way as r — s- oo. This result and the method of proof 
will be used in the analysis of solutions of the afore mentioned boundary layer problem. 
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>, 

t^ ■ 1 Introduction. 



The analysis of thin films coating surfaces is a basic problem in fiuid mechanics and has many 
applications in industrial processes. In particular flows on surfaces with topography had been 
extensively studied (see, for instance, [S], [H], [13]: HH, [E], [IS], [10]) ■ Usually, the presence of 
regions where the liquid accumulates plays an important role in the overall dynamics of these fiows. 
In many situations, the arrival of a thin layer of fiuid to a region where the liquid accumulates 
or reduces its velocity, results in the formation of oscillations on the fiuid interface (see e.g. [2], 
[1], [2], [21] and references therein). On the other hand, it was claimed in [18] that the pressure 
difference between phases in unsaturated porous media fiow, when complete wetting is assumed. 
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3 . is mostly due to the regions where the fiuid accumulates. 

Our goal is to derive precise mathematical results concerning the liquid accumulation regions 
mentioned above. To this end, we will study a particular type of fiows. Namely, we are interested 
in the motion of thin layers of liquid moving on an smooth and, for the moment, arbitrary surface. 
Specifically, our goal is to understand the dynamics of such layers of liquid in the case in which 
the forces induced by the curvature of the substrate are comparable to gravitational forces. For 
this particular rescaling of the force terms, the liquid tends to accumulate in the regions where 
the gravitational and surface tension forces balance. Outside such regions it is possible to simplify 
the description of the fiuid motion using a thin-film approximation. However, this approximation 
breaks down in the regions where the difference of the forces induced by gravity and surface tension 
vanish, these are the regions where the fiuid tends to accumulate. As a consequence, a boundary 
layer is required to describe the fiow there. 

In this paper, we will restrict our attention to stationary fiows in two space dimensions. In 
this case, the thin-film approximation for steady fiows reduces, to leading order, to 

§-^{Qis)h')=0 (1.1) 
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where h is the (non-dimensional) height of the fluid above the substrate and Q{s) is a function 
that depends on the geometry of the substrate (parametrised by its arc-length s) and its relative 
orientation with respect to the direction of gravity. 

As long as the function Q{s) has a constant sign it is possible to obtain a stationary solution 
of (|l.ll) for any given flux describing the flow. However, if the function Q{s) changes sign it is not 
possible to obtain steady flows satisfying globally the assumptions required for the validity of the 
thin-film approximation. We remark, that it is possible to assume in this case, using a convenient 
parametrisation of the curve where the fluid lies, that Q{s) > 0. We shall assume this condition 
in most of the regions where the fiow takes place, except in a very small region where the liquid 
will tend to accumulate. In such region the difference of curvature between the substrate and the 
interface, that has been neglected in the derivation of (|l.ip . becomes relevant. Nevertheless, the 
thin-film approximation remains valid in all of the flow considered here. 

We obtain two main results. First, we derive, using asymptotic approximations, an equation 
describing the height of the liquid in the (small) accumulation regions where Q{s) < 0. This 
equation reads, in suitable non-dimensional variables, 

'^'^ e + a)H^ = l, (1.2) 



where H is a rescaled fluid height and ^ is a rescaled arc-length. The number a G M is a parameter 
that provides information about the size of the region where Q{s) < 0. Equation p.2p must be 
solved subject to the corresponding matching conditions, namely, |^|3 as |^| — >■ oo. The existence 
of such solutions is not a priori clear, and their existence will be proved in [9]. The solutions 
obtained there are the natural candidates to describe the boundary layer arising at the stagnation 
or accumulation points, i.e. the points sq where Q{so) = 0. In the analysis of p. 21) it will be 
essential to have a precise description of the solutions of the following auxiliary equation: 

(0.i)*..i. ,o) 

Observe that (|1.2p is a non-autonomous third order ODE. Equation (|1.3p results when we introduce 
new variables that transform (|1.2p into a fourth order autonomous system, and that make the 
matching conditions also autonomous. Namely, in the new independent variable r, (jl.2p reduces 
to (|1.3p as \t\ — > cx). That means that the limits r — > — oo and r — > cx) define two invariant 
subspaces for (|1.2p where the flow (of the dynamical system) is given by (|1.3p . In these variables, 
the solutions of (|1.2p that satisfy the matching conditions correspond to a heteroclinic connection 
(the critical points connected correspond to the matching conditions). The existence proof of 
these heteroclinic connections is done by means of a shooting argument and requires detailed 
information about the flow given by (|1.3p . that is the flow near \t\ — >■ cx). 

Our second result in this paper thus concerns the rigorous analysis of (|1.3p . As we shall see 
later, one of the main features of (|1.3p is the existence of solutions that oscillate with increasing 
amplitude as t — >• cx). We thus investigate this oscillatory behaviour rigorously in this paper. 

Equation (|1.3p has an independent physical interest. It arises in the description of a thin-film 
free surface flow on a vertical wall that is moving downwards with constant velocity (cf. [23]). 
The existence of oscillatory solutions with increasing amplitude was observed in [23j using formal 
asymptotic. These oscillations are related to the possibility of air entrainment in these flows. Due 
to the fast increase of the amplitude of oscillations, the validity of the thin-film approximation 
breaks down and the corresponding solutions need to be matched into other flows. In fact, in [23j 
the solutions are matched into a stationary pool of fluid. 

We observe that the onset of oscillatory interfaces with increasing amplitude seems to be a 
general phenomena of thin-film flows arriving at a region where the flow is slower. For example, 
a similar phenomenon is observed in [7] , where the motion of air bubbles in a fluid that is flowing 
inside a tube is analysed. When the tube is placed horizontally and its radius is supposed to 
be small enough as to neglect gravity effects, equation (|1.3p . but with a power 2 non-linearity, is 
derived to describe the flow of fluid separating the bubble surface from the tube wall. This flow 



is then matched into a minimal surface describing the front tip of the bubble. We point out that 
oscillatory behaviour for the equations obtained in [7] and in |23| have been numerically reproduced 
in |21| . Another example of such oscillatory behaviour of thin- film flows in the presence of a slower 
flow ahead is provided by thin- film flows on inclined planes (cf. [1], [5]). 
We notice that general third order nonlinear ODEs of the form 

^ = /(^), (1.4) 

arise in many other problems of fluid mechanics and in the present context of thin-film approxima- 
tions in particular (cf. [10], [11], [15], [16], [17], [22]). Some of them have been studied rigorously. 
Typically, the type of problem that has been considered in the rigorous mathematical literature is 
the possibility of having solutions that vanish at a finite value of r, as r increases. For instance, in 
[3] a necessary condition on the non-linearity /, such that $ vanishes in the forward direction is 
derived. In the particular cases /($) = 1 — 1/$" with n = 2 and n = 3, this necessary condition 
is not satisfied. This implies that the solutions of ()1.4p are globally defined and remain always 
positive for arbitrarily large values of r. Qualitative properties of the solutions of the equation 
(|1.4|) with a non-linearity of the form 

/($) ^^— + ± 

have been studied in [6]. Some specific non-linearities relevant to applications in fluid mechanics 
are described in [3], [6] and |21j . In this last paper, numerical solutions showing oscillatory 
behaviours for some class of non-linearities are obtained. However, to our knowledge, no rigorous 
result describing the oscillatory behaviour for the solutions of (|1.4p have been given before. 

The paper is organised as follows. In Section [2] we derive (|1.2p . We perform a thin- film 
approximation on the corresponding Stokes flows, that in two dimensions, is written in curvilinear 
coordinates. We derive the leading order problem and analyse the significance of the driving term 
Q{s). As it was pointed out earlier, this coefficient is completely determined by the geometry of the 
substrate. In order to justify our assumption on Q(s), we also prove in this section the existence 
of curves for which Q{s) has the asymptotic behaviour Q{s) ~ As"^ as s — ^ for some A € M. 
We end the section by deriving the boundary layer problem near a point where Q{s) vanishes, 
that is, we derive equation ()1.2|) and the corresponding matching conditions. The expressions in 
curvilinear coordinates of the linear operators needed to derive the next order term in the thin-film 
approximation are given in Appendix 1. 

In Section [3] we analyse (|1.3p . We start in Section [3Tl with the local behaviour near the crit- 
ical point and the behaviour of solutions on the (one-dimensional) stable manifold, and the ones 
outside this manifold. For the former ones we prove that either $ — >■ oo as r — >■ — cx) or $ — >■ O"*" 
as T — > (r*)+ for a finite r*. For the later ones, we prove rigorously that there exist a two dimen- 
sional family of positive solutions of (jl.3p . defined for — oo < r < oo, satisfying limT-^_oo ^{t) = 1, 
limsup^_j.o^ (f>(r) — oo and liniinf^_j.oo *&(''') = 0. These solutions thus oscillate with increasing 
amplitude, a fact observed and described in [23]. In Section [3?2l we make precise and prove rigor- 
ously the asymptotic formulae obtained in [23] describing the rate of increase of the oscillations. 
Moreover, we prove that this precise type of oscillatory behaviour is an attractor for the solutions 
of (|1.3p . in the sense that every solution of this equation that does not approach $ = 1 as t — > oo, 
eventually oscillates in that precise way. 

Finally, we remark that the oscillatory behaviour of solutions of (|1.3p will be shown to arise due 
to the transition between ranges of r in which p.3p can be approximated either by the equation 



when $ is large, or by the equation 






1 =0, (1.5) 



$^ = 1, (1.6) 



when $ is very small. Equation (|1.5p can readily be integrated. The analysis of (jl.6|) is done 
in Appendix 2, where the relevant information about the solutions is gathered in a theorem. 
Equation (jl.6|) is invariant under the transformation (r, $) — >■ (Z^r, /$) for any Z G R and can 
thus be analysed by means of a related transformation that reduces the study of this third order 
equation to the study of a planar autonomous system of ODEs. 

2 Derivation of the model 

We first derive the equations describing the motion of a fluid in a stationary regime flowing 
down a substrate. Since we restrict ourselves to two-dimensional flows, the liquid phase fills a 
two-dimensional domain fl. We assume that the fluid wets completely the surface. Then, the 
boundary of this domain can be decomposed in two disjoint pieces dfls and dflf. The portion of 
the boundary dQs denotes the intersection of the liquid with the fixed substrate. On the other 
hand, di^f is the free interface separating the liquid and the gas phase. 

We let he denote the characteristic height of the liquid over the surface, and Rg is the typical 
radius of curvature of dflg. We henceforth assume that j^ is small in most part of the flow. This 
allows us to use a thin-film approximation. In particular, this means that the normal component 
of the velocity is very slow compared the to tangential one. 

We deal with very viscous flows, thus we assume that the flow takes place at a very low 
Reynolds number (Re ^ 1) and, therefore, that inertia effects arc negligible compared to viscous 
forces. Taking into account our initial assumption, the Reynolds number is Re = h'^vt/{Rsv)- As 
it is usual in the study of viscous flows under the thin-film approximation, the order of magnitude 
of the tangential velocity along the surface is 

..«^, (2.1) 



where v is the kinematic viscosity. The total flux of fluid per unit length along a surface transversal 

t: 

gh'^ he J 



to the substrate is then of order J « S—^. With this, we observe that in particular 
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i.e. the main assumption on the ratio of the characteristics lengths allows rather jast very viscous 
flows (cf. [J). 

Under these assumptions we can describe the flow using the free surface Stokes equations. In 
Cartesian coordinates (x, y) in ID we take the gravity vector pointing downwards and so g = ^9^y 
where e^ = (0, 1)^ and g is the gravity constant. 

We non-dimensionalise the problem with the length i?^, as follows: 

g = -g^y , X = i?,x* , u = ^u* , V = pgRsP* , K = ^K* , t = ^ t* , (2.2) 

1/ Rs g Rs 

where u is the velocity field, p the fluid pressure, t is time and K is the curvature of the free 
surface (951/). In these non-dimensional variables, where we drop the * for simplicity of notation, 
the Stokes flow reads 

Au - Vp - Gj, = , divu = , in D,, (2.3) 

and is supplemented with kinematic and surface tension boundary conditions on the free surface: 

VN = u-N , (Vu + (Vu)^)N -pN = —KN , on dQf , (2.4) 

B 

and with the no-slip boundary condition on the substrate: 

u = on dfls- (2.5) 



Wc notice that the order of magnitude of the velocity is given by (j2.ip . Since in problems where 
the thin-film approximation is valid the main component of the velocity is the tangential one, using 
(|2.2p . it then follows that the order of magnitude of the (non-dimensional) tangential component of 

u is ( -^ j <C 1. The other dimensionless parameter appearing in (|2.3p - (|2.4p is the Bond number: 

^ ^ W_R? _ P9 Rs 



a/Rs ' 

that measures the relative size of gravity and surface tension or, more precisely, the ratio between 
the hydrostatic pressure due to changes of height of order Rs and the changes of pressure induced 
by the surface tension on surfaces with curvature radius Rs- 

In this paper we assume that B is of order one. In particular, since |u| ^ 1, the second 
equation in (|2.4[) implies that p ^ ^ K, i.e. curvature and gravity effects are of the same order of 
magnitude; that is the characteristic feature of the limit considered in this paper. 

2.1 Curvilinear coordinates 

We now proceed to derive the equations describing the height of the liquid using the thin-film 
approximation. To this end we reformulate (|2.3p - (j2.4p using a more convenient set of curvilinear 
coordinates. 

Let the substrate be described by a regular curve parametrised by its arc- length s: 

dns = {a{s) = {x{s),y{s)) : s G R} . 

We let n denote the outer normal vector of dfts (thus when the substrate is horizontal and s 
increases, if {x{s), y{s)) moves towards the right then n points upwards), i.e. n(s) = {—y'{s), x'{s)). 
Then the curvature of dV,s is k{s) = x'{s)y"{s) — x"{s)y'{s). We shall also denote the tangent 
vector to dV,s by t(s) = {x' (s) , y' (s)) (observe that then det (t(s),n(s)) = -1-1). We introduce the 
coordinates (s, d) such that a point in space is determined relative to the substrate by its distance 
to it: 

X = {x{s, d), y{s, d)) = a{s) + dn{s) . 

The corresponding orthonormal basis is then {ei = t, 62 = n}, which only depends on s. We 
shall let the vector components in curvilinear coordinates have indexes 1 (component tangential 
to the substrate) and 2 (component normal to the substrate). 

The free boundary is described by an unknown function measuring the distance to the substrate 
d — h(s,t), hence it is parametrised by 

X/(s) = a{s) + h{s)n{s) , 

its tangent and normal vectors are given by 

, , . dh 9/)- , , , , 

{l-kh)t + —-n -—-t + {l-kh)n 

T = Tit + Tan = ^ — , N == iVit + iVan = Q^ 
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The curvature of the free boundary is 
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r. _ os^ OS'' \osJ ds ds 
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A transformation of the relevant differential operators to the new coordinates is given in e.g. 
[T]. Using this, problem (|2.3p - (|2.5p in the new coordinates becomes 

(2.7) 
L — Ku \ us J ua 

and 
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(2.9) 



with 



((Vu + Vu'^)N) •N-p= — A' on d^h 
((Vu + Vu^)N) •T = on d^h 

where the full expressions of the tangential and normal components of (Vv + Vv^)N are given in 
Appendix 1. Finally, the second equation in (|2.4p becomes 



^^ ui dh 

m='''-T^kdTs °" '^^'^ ('-^'^ 

We also observe that using this coordinate system, the flux J per unit length along any surface 
orthogonal to the substrate, can be written in the original dimensional set of variables (u, x) (see 
(1221)) as: 

J= ^^-^ / uidy. (2.11) 

'^ Jo 

Obviously, J^ ui dy is just the non-dimensional flux. 

2.2 Thin-film approximation 

We will assume that in most of the fluid we have: 

h, 



«1. (2.12) 



Combining (|2.ip and (j2.2p it then follows that, in the non-dimensional variables, ui « e^ ^ 1, 

therefore, also that L uidy w e^ ^ 1 (and this in particular means, in view of (|2.1ip . that 
Ji^ ~ ^3 ^^ 1 ^ ^jr„ i- 



-^ ~ e^ <^ 1). We the rescale the variables as follows: 



rd* , h = eh* , t = et* , ui ~ e u^ , u^ = e U2 , (2.13) 



where the last one results, to leading order as e ^ 0, from the rescaling of (|2.7p . With these 
changes, (J2.7p - (|2.9I) reduce, to leading order as e -^ 0, to 

M + M-o ^'"1 ^P c -t-O -^-0 (2 14) 



with boundary conditions 

dd* "' " B 



^^0,-P=1k, on d*=h* (2.15) 



where, from 



K^k + e(^+k^hA+0{e^). (2.16) 

The no-shp condition p.Sp stays as is, and p.lOp becomes 

to leading order. 

These equations are now combined into a single one for h. Integrating the second equation in 
(I^Ti)) one gets 

I / dp 



and, from the first equation in (|2.14p . one gets "2 = ^/0 IT^ ^V^ ^^^ 



Equations similar to (|2.15p . (|2.16p and (|2.18p . where the main driving terms are the gravity 
and the curvature of the substrate, have been obtained, in a slightly different context, in [T7] and 
[TB] . This model can be obtained also as a particular case of the ones considered in [TH] for specific 
choices of the parameters. 

We shall write the leading order of (|2.18p as 



f -sJil*) ("■)') ^°. -- «<"-(5S 



1J77 +o-(Q(«) (/**)') =0' ^ith Q(^)^= ( o— -ej^-t). (2.19) 



This is a non-linear non-homogeneous (in s) hyperbolic equation. The term Q (s) measures the 
tendency of the fluid to move in the direction of the tangent vector t. It is important to remark 
that the function Q (s) is completely determined by the geometry of the curve and the gravitational 
field. Moreover, all the geometrical features of the curve are reduced in this approximation to the 
function Q{s). 

We are interested in the study of this system in the steady state regime. Which, to the current 
order of approximation, is simply 

^Q{s){h*f = l. (2.20) 

We observe that (|2.20p can have solutions with Q{s) < 0, h*{s) < 0. The meaning of negative 
values of h* is just that the fluid is placed in the direction of the vector — n(s), i.e. if a point 
moves along the curve dClf in the direction of increasing s, positive values of h* , indicate that the 
fluid would be seen by the moving particle to the left , and negative values of h* mean that the 
particle would see the fluid to its right. We can always assume that Q{s) > by reversing the 
direction of the parametrisation. 

We consider the case in which Q{s) vanishes for some values of s. At these points the ap- 
proximation (|2.20p would break down, because the height of the fluid is expected to increase 
unboundedly; the next order terms (in particular, in (j2.16p ) will become important near these 
points. We will not consider in this paper curves for which Q{s) takes negative values over some 
intervals. 



2.3 Geometrical problem of the substrate 

We are interested in studying flows for which the approximation (j2.20p breaks down; that is, if 
Q(so) = for some value of sq, where, for definiteness, Q{s) > for s y^ sq in a neighbourhood of 
So- It is then relevant to show that there exist surfaces where these conditions are satisfied. We 
prove in this Section that this is indeed the case. 

We begin formulating the problem satisfied by curves if Q{s) is a given function. We recall 
that, since a{s) = {x{s),y{s)) is the arc-length parametrisation of dils, {x'{s))'^ + (j/'(s))^ = 1. 
Let us denote by 6{s) the angle between the tangent vector t and the horizontal axis. Then: 

^=cos(0(s)), ^=sin(0(s)). (2.21) 

as as 

Observe that (restricting to values oi 9 G [— 37r/2, tt/2]) when 9 G (— 37r/2, — 7r/2) the fluid is below 
the substrate, and, when 9 G (— 7r/2,7r/2), then fluid is on top of the substrate. For 9 = —tt/2 
the surface is vertical and the fluid is to the right of it. Similarly, for 9 = tt/2 (or —3tt/2), the 
substrate is vertical and the fluid is to the left of it. In terms of 9, the curvature is given by fc = ^ 
and By ■ t = sin {9{s)). We can now write (|2.19p as 

l^-sm9{s)=Qis). (2.22) 

This is the differential equation that describes the geometry of the substrate in terms of the 
function Q{s). Solving (|2.22p the curve is recovered by integrating (|2.2ip . Observe that (|2.22p is 
the equation of a forced nonlinear pendulum. 

There is a class of semi-explicit solutions of ()2.22p that can be obtained assuming that Q (s) = 
Qo G IR- The simplest case corresponds to IQol < 1, and includes the pathological case Qq = 0, 
for which the thin- film approximation breaks down everywhere on the substrate (see (|2.20p '). For 
such Qo's there are constant solutions of (|2.22p given by 

9{s) = 9o , sm{9a) = -Qo ■ (2.23) 

These solutions correspond to an inclined plane with constant slope 6*0. We can always assume 
without loss of generality that Qq > 0. We will assume also that 6*0 S [— 7r,0]. Due to our choice 
of signs, the fluid is above the plane {9 (s) = ^o} if ^o G (~f : O] and the fluid is below the same 
plane if 9q e [— tt, — f ). In the critical case 6*0 = —§ the plane is vertical and the fluid is to the 
right of it. Also, for any < Qo < 1 there are two roots of the second equation of (|2.23p with 
^0 G [— TT, 0]. Let us them be denoted by 9i and 9^, where 

-Tr<9l<-^<9i<0, 

and they satisfy 91 + 9i = — tt. 

We can also consider perturbations of these constant solutions. The point (6*, ^) = {91, 0) is a 
centre in the phase portrait of the equation (|2.22p . Therefore, there is a one-parameter family of 
periodic solutions of (|2.22p for which 9 (s) oscillates around 91 in the usual way for the nonlinear 
pendulum. We can then obtain curves dfig by means of (|2.2ip . The resulting curve does not have 
self- intersections if the amplitude of the oscillations of {9 (s) — 91) is small, as it can be seen using 
a continuity argument. It is interesting to notice that for all these surfaces the height of the liquid 
remains constant (sec (j2.20p ). even if they are not planar. The fiuid is below them, but we would 
assume that the fluid is stable enough to allow such flows. 

We remark that the quantity H {j^,9^ = -j^ (^) + cos{9) + sin(6'o)0 is conserved under 
the flow (|2.22p . Using this, it can be seen that there is a homoclinic orbit connecting the point 
(^' ^) ~ (^1' '-') ^° itself. In view of (|2.20p . this provides another example of a surface on which 
the height of the fluid remains constant. A perturbative argument shows that the corresponding 
surface obtained by means of (|2.2ip approaches asymptotically, as s — ?► ±oo, the plane {9{s) = 9i} 



and it docs not have self- intersections if 9q is close to — ■^. However, since this is not the main 
goal of this paper, we will not continue this discussion here. 

A set of interesting solutions of (|2.22p are the ones associated to \Qo\ > 1. For these solutions 

\0{s)\ increases very fast as s -> oo. Their asymptotics are 9{s) ^ 'y^ as s — ?> ±cxj. The 
corresponding surface obtained by means of (j2.2ip approaches asymptotically a finite point in each 
of the limits s — >■ cx) and s — >■ — cx). Therefore, these surfaces cannot be extended to unbounded 
domains and the approximations yielding to p.20p must breakdown at some point. However, these 
curves give interesting examples of substrates where the driving coefBcient Q{s) is constant and 
larger than the maximum value allowed by the gravitational force, due to curvature effects. 

After these preliminary observations about (|2.22p , in the following theorem wc construct func- 
tions Q{s) that vanish quadratically at sq ~ 0. 

Theorem 1 Let < Qo < 1- There exist ujq > small, such that for any uj < wo? there exist 
curves dfls € C°°(R) without self-intersections that approach asymptotically, as s ^ ±oo, the line 
{9{s) = 9i}, where 9i e (^f lO) and sm.{9i) = —Qo, and such that the function Q{s) defined by 

means of WJ^) satisfies Q{s) >Q for s^Q, Q{Q) ^ uj and < lim^^o ^^"jr" < oo. 

Proof. Wc construct the function Q(-) as follows. We take Q{s) ^ Qo (l — £, (£;)) for s < cj, where 
C e C°°{R), Ciz) > for z e (0,i), Ciz) < for z e (i,l), ^(z) = for z"e (-oo, 0] U [1, c3o), 

^" (^) ^ *-* '^^'^ /o ^i'^)'^^ — 1- The solution 9^ of the equation (|2.22p can then be approximated 
for s < w, if w is small enough by means of the solution of the distributional equation: 

1 d'^0 
— — ^ - sin 6l(s) = Qo - Qo^S{s) . 

An alternative way of proving this can be obtained by rescaling s with u and studying the resulting 
regular perturbation problem using Gronwall. We then obtain: 

lim'^"(-)'^^'=0, u^l^^(-)-^Qo-l=o 

w— >0 bJ w— >0 U) 

We now construct 9i^{s), but not Q{s), for s > w as any function that extends the obtained 
9^{s) in the region s < w as a function in C°°(R) with \9^{s) - 9i\ < |, \9i{s)\ < 2BQolo, 
|^"(s)| < 2Buj for s > id. We impose also that Mms^co&uis) = ^i- This extension is possible 
because |6'"(s)| ^ w as w ^- 0. We then define Q{s) by means of (|2.22p for s > w and construct 
the curve dflg by means of (|2.2ip . Since \9ui{s) — 9i\ < j ior s £ R it follows that the resulting 
curves do not self-intersect for ui small enough. The rest of the properties stated in the Theorem 
are straightforward. ■ 

We recall that our sign criteria implies that the fluid lies above the curve dftg obtained in 
Theorem [2 

2.4 Analysis of water accumulation regions in the limit e — t- 

We now study the form of the stationary solutions of the thin-film approximation in the regions 
where the leading order approximation (|2.20p breaks down. We will consider only the ease in which 
Q"{sq) > at the points s = so where Q{s) vanishes. The existence of such non self-intersecting 
curves dils with Q e C°° (M) satisfying (j2.26p below has been proved in Theorem [TJ Assuming 
(|2.12p and performing the scalings (|2.13p into the Stokes problem (|2.8p - (|2.10p . one can derive the 
stationary thin- film approximation to order e, namely, dropping the *'s for simplicity of notation. 



with 



p^, \ d ^2 7, ^, , 5 1 dk 3d. 

'^^'^ = Ts' -s'^^^'^^sBd-s-sa-s^''^ 



and Q{s) as in (|2.19p . Assuming now that the flux is the same as the one obtained in the leading 
order approximation (|2.20p . wc obtain the stationary equation: 



i^^h r)h 8h ~ \ 

Qis) + e^ + efc' ^ - eiey ■ n)^ + eQis)h U^ = 1 . (2.25) 

os-^ as OS J 



We take a somewhat general form for Q(s), namely, we assume that near a stagnation point 
placed at s = we have: 

Q{s)=As^ + be^+o{s'^+e-r?), A>0, beR. (2.26) 

The exponent j= results from the fact that with this particular rescaling the effects of this term 
will turn out to be of order one in the boundary layer. We will also assume that Q (s) > for 

|s|>ei?. 

The natural rescaling of variables required to study the singular perturbation problem (J2.25P 
is: 

Finally, (|2.25p then becomes (|1.2p . to leading order, where 

f 2 

a 



V3Ai 
In order to determine the boundary conditions that must be imposed to the solutions of (|1.2p 



we must study the asymptotics of the solutions of ()2.20p as s ^ 0. Using p.26|) . we obtain from 

h{s) ^ (-^\ ' as e^ « \s\ « 1 , (2.28) 

and using (|2.27p . we find 

JJ^J_ as l<|e|<£"*, (2.29) 

which provides the matching conditions for the solutions of (jl.2p . 

According to the picture emerging from the results described here it would follow that the 
stationary solutions of (|2.24p could be approximated by means of the solutions of (|2.20p (i.e. by 
(|2.28p ) except near the stagnation points or where Q{s) is very small. At such points, a boundary 
layer arises and it is described by means of the rescaling (|2.27l) and the function H that solves 
(fO]) with matching conditions ()2.29p . 



3 The solutions of ( ITTSi ) 



In this section we analyse (|1.3p rigorously. This analysis will be used in the proof the existence of 
solutions of p.2p subject to (|2.29p that is done in [3]. 

3.1 Global existence, stable and unstable manifolds 

In this Section we describe in detail the solutions of (|1.3p . We write (|1.3p in the equivalent system 

We observe that there is a unique critical point for (|1.3p . namely P, = (<i>,VF, ^) = (1,0,0), 
and is hyperbolic. The stable manifold of (|1.3p at the point Pg is tangent to the vector: vi = 
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(3 3,-3 ^,1) , and the corresponding eigenvalue is Ai := —33. Its unstable manifold is two- 
dimensional and tangent at the point Ps to the plane spanned by the vectors V2 '■= (— g36,g33,lj 

and f3:=(i35,i3'5,0) . The corresponding eigenvalues of the linearised problem being complex 

conjugates, namely, A2 := 33 ( 1 + i 33 j /2 and A3 := 3^ ( 1 — i 33 j /2. 

After these preliminary observations, we next prove that the solutions of p. 31) do not develop 
singularities for increasing r. 

Lemma 2 (Forward Global Existence) Suppose that $(to) > 0, W{to), ^(tq) are arbitrary. 
Then, the solution of LS. 1]) with these initial data is defined and $(t) > for any t > tq. 

Proof. The only possibility of losing global existence is when $(t) approaches zero at a finite 
value of r, because when $(t) > (5 > a simple estimate yields, integrating the equation, that 
($, W^,*) < (p(z),p'(r),p"(T)) where p(t) = $(ro) + ^{l/S^ - 1)(t - Tof + vl/(ro)(r - to)V2 + 
W^(ro)(To), thus the solution cannot become unbounded at a finite value of r. 
First, we show that if 

lim inf ^ $(t) =: (3.2) 

then 

lim $(r) =0. (3.3) 

We prove p.3p by contradiction. First, we observe that as long as the solution is defined 

d^ 1 

holds and then ^(t) > ^(0) — r, therefore 

*(t) > -Ct* for < T < T*and some C > . (3.4) 

Let us now assume that (|3.2p and that linir-s-r* ^{t) does not exist. Then, for some small enough 
Eq > 0, there exists an increasing sequence {t„} such that lim„_j.oo Tn ~ t* , $(t„) ~ 2£o and 
W{Tn) > for aU n e N. Lrtcgrating (|3.ip gives 

$(r) = $(r„) + W{t„){t - r„) + / / *(s)dsd77 , for t > Tr, , 

hence, using p.4|) . we obtain that, as long as ^(r) > Eq, 

$ (r) > $ (t„) + M^ (t„) (t - T„) - / / /^T* ds dr/ 

K ^ 2 

> 2£o - yT* (t - r„) , for r* > r > t„ . 



It then follows that $(r) > Eq for t„ < t < r„ + w 3°^x . On the other hand, due to p.2p . there 
exist {tVi} such that lim„_j.oo -r„ = r* and $ (f„) < ^. We can assume that f„ > r„. Since for 
n large enough we have f„ € (T„,r„ + J^f^]i we obtain a contradiction, and therefore 



must hold if p.2p does. 

We now prove that (|3.3p cannot happen. In order to make this precise we employ the change 
of variables $ = (r* - T)^i^, W = (t* - t)^W , ^ = (t* - T)~i^ with s = - ln(T* - t), which 
gives the system 

rf$ ,7, 3- dW - 1- d* 5-1 _9, ,^^, 

ds 4 ds 4 ds 4 $3 
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Wc observe that with this change of variables, when r — > (t*)~ then s — > +00 and the inhonioge- 
ncous term in (|3.5p becomes very smaU. In the absence of this term, the remaining autonomous 
system has no critical points, suggesting that all trajectories are unbounded. 

Let us assume that (j3.3p holds, then, in terms of the new variables, this is equivalent to 

lim e~3-'$^o. (3.6) 



First we assume that lims_j.oo 'J' 7^ 0. Then, multiplying the third equation in p.Sp by e^'^ we 
obtain: 

thus for large enough s, (|3.6p implies that there exists a constant C > such that 
Setting s > sq with sq large, integration of p.7p gives 

Jso 

hence, taking for example ^0 = e^^''^(so) for any such sq, we obtain that for s > sq there exists 
a constant (of undetermined sign) such that ^(s) > e^^'^^'o- Proceeding in a similar manner, the 
second equation in p. 51) yields W > e~'i''Wo for large s, and some constant Wq. Using this now 
for the first equation in (|3.5|) gives, for large sq and s > sq, 

e-^'^ > e-i'°^{so) + e-'^Wo - e'^Wo . (3.8) 



Since $(s) does not tend to as s ^ 00, one can always choose sq large enough such that 
$(so) + e~4*''Wo > S > 0, and then (|3.8p implies that lim^-j-oo e"*"!" > 0, a contradiction. 

If lims_j.oo $ = 0, we argue in a similar way directly on the equations (|3.5p . The third and the 
second equations give that there exist a large positive constant C such that ^, W^ > C for s large 
enough. Then the first equation in p.Sp and the fact that i'(s) > or all s imply that ^$ > C, 
a contradiction. ■ 

The next lemma deals with the solutions contained on the stable manifold of the critical point. 

Lemma 3 (Behaviour of solutions on the stable manifold) // ($, W, ^) ^- P, as r ^- 00 

and ($, W, *) 7^ Ps, then either 

$ -^ +00 as r — > -00 (3.9) 

or there exists a finite t* such that 

$(t)^0 as t^{t*)+ with $(r) > for all t > t* . (3.10) 

Proof. We first observe that integrating (|1.3p three times we obtain 



''""^^- '*°°' ' iVn, ^^rrfT^i^-lMn*. (3.11) 



dr' J, \mn)f j '• dT J, ;„ V(*(ri))' 

(•+00 /.+00 /.+00 / 1 \ 

*W-1 = / I I (^^--—^dr,dr,dr,. (3.12) 

Due to the fact that the trajectory under consideration is in the stable manifold of Pg the 
corresponding orbit of ($, W, \E') is tangent to Ps in the direction vi. Suppose that this happens 
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in such a way that there exists a f large enough such that $(t) > 1 for all r e (f, +00). Let f be 
f = infJT : $(cr) > 1 V cr S (r, +00)}. 

Our goal is to prove that f = — cx). Suppose that f > —00. Then, the continuity of $ implies 
that <E> (f) = 1. Using (|3.12p with t = f and the definition of f we obtain a contradiction, whence 
f = — OD. Then p.ll[) implies J-^' < C < for any t sufficiently small, and therefore p.9p 
follows. 

Suppose now that the orbit ($, W^, ^) is tangent to Pg in the direction —vi, then there exists 
f such that $(r) < 1 for all r e (f,+oo). Now ((3lT|) implies that ^ > for all r e [f,+oo). 
Thus, trivially, $(t) < 1 for all r G (t*,+03), where r* is defined by t* = inf{r : $(cr) > 
V (T G (r, +00)}. Let us prove that t* > —00, so that in particular, by continuity, $(r*) ~ and 
$(t) e (0, 1) for all T G (r*, +00). If, to the contrary, r* = -00 then, since ^ > and ^ < 0, 
it follows that$(r) < $(ri) + $'(Ti)(r - ri) for any ti G R, hence r* > -00. ■ 

The solutions on the stable manifold are all those solutions with $(r) ^ 1 as t ^- 00: 

Lemma 4 // $(t) is a solution of il.c^) with liuir^ao $(t) = 1, then ($, P^, *) ^ Ps . 
Proof. We define P(r;ri) as: 

1 



VI rV2 



P(r;ri) = <I>(r)- / / / 



i> (%)' 



1 



drjs dri2 dr]i 



(3.13) 



Then, using lim^_j.oo ^{t) = 1 as well as (JSHJ it follows that for any ti, P{t; ti) is a second order 
polynomial given by 



P{t- n) = $(ri) + — f^(r - n) + —^{t - T,f 



dr 
Taking ri — > 00 and using (j3.13p . it is clear that 



dr2 



(3.14) 



lim P(ri; Ti) = lim P(ti + 1; n) = lim P(ti + 2; ti) = 1 

ri— foo ri— >-oo ri— >-oo 



Therefore p.l4p implies: 

"d$(Ti) d2$(^^) 



lim 



dr 



dr2 



0, 



lim 



d$(Tl) _ 2.^'*(^l) 



dr 



dr^ 



0, 



and hence ($, W^, ^) ^^ P^ as r ^ 00. ■ 

Wc now give the basic behaviour of the solutions that are not in the stable manifold: 

Lemma 5 Suppose that ^ is a solution of il.3]} such that $(t) does not converge to 1 as t —> 00. 
Then: 

lim sup ( —2 +'!>)= +00 . (3.15) 



Moreover, 



2$2 



lim sup $(t) = +00 and lim inf $(t) = 0. 



(3.16) 



dOl) by ^ we obtain: 



Proof. There is a monotonicity property associated to the solutions of p.ip . Indeed, multiplying 

dr 

d {d^<^d<^ ^ 1 ^^ _ (d^'^^ 
^ ^ ^ Ih^ 



hence 



dr ydr'^ dr 2$2 
*2 > , E 



>0 



dE 
d7 



*VK 



1 
2$2 



$. 



(3.17) 
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The energy estimate p.l7p yields estimates for ^^ for the trajeetories eontained in the unstable 
manifold. More precisely we have: 



{^{s)fds = ^(oo) - EiTo) 



for any tq E M. Notice that E{oo) could be infinite. Actually, our goal is to show that E{oo) = oo 
unless $ — > 1 as T — > oo. 

Suppose that E{oo) < oo, then 



{'<S/{s)) ds < oo. 



(3.18) 



We claim that if this is the case then limT-_j.oo ^(t) = 0. Indeed, the third equation in (|3.ip yields: 



— — > -1 for ah r G 

dr 



(3.19) 



Suppose that there exists a sequence {t„} such that lim„_^oo t„ = oo and satisfying ^(t„) > Sq > 0. 
Then (|3.19p implies ^(r) > Eq — {t — Tn) for all t > Tn, hence 



T +^ 



(^(T))^dT> ^ 



and this gives a contradiction with (|3.18p . Similarly, suppose that there exists a sequence {t„} 
such that lim„_j.oo t„ = oo and satisfying \I'(t„) < — Eq < 0. Then, using again p.l9p we obtain 

^(t) < -£o + (t„ - r), r < t„, hence 






(*(r))'dT>^, 



that also yields a contradiction, therefore 



lim *(r) = 0. 



(3.20) 



The first two equations of (jSTTj) imply: 

sup \W (s) - W {t)\ ^ as T ^> oo 

sG[t,t+1] 

sup |<I>(s) - $(t) - M^(t)(s - r)| ^ as r ^ oo , 

sG[r,r + l] 

On the other hand, integration of the last equation of p.ip implies 

nr+l 



*(r + l)-*(T) = / 
which together with p.20p gives that 



1 



($(.))^ 



- 1 



ds . 



lim 

r— )-oo 



r+1 



ms)y 



ds^O. 



(3.21) 

(3.22) 



(3.23) 



Therefore (infgg[^^_|_i] (f>(s)) < 2 for r sufficiently large, since otherwise there would be a 
contradiction. It then follows that lim-r- i-oo M^(''') = 0. Indeed, otherwise (|3.2ip . p.22p would 
imply that $ would become negative for large values of t. It then follows from (|3.22p that: 

sup |<I>(s) - $(t)| ^ as T -> oo 

se[T,r + l] 
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and (|3.23p then yields liniT-^oo ^(t) = 1 against the hypothesis of the lemma. The contradiction 
yields £'(00) = 00. 

We now prove p.lSp . Suppose that 

1 
2$2 



lim sup ( — ^ + $ ) < 00 . (3.24) 



Since E{oo) = 00, it follows from the definition of E in (I3.17P that 

lim [*(t) W{t)] = 00 . (3.25) 

r— foo 

Suppose that limj_j.oo |^(7j)| = 00 for some sequence {tj} with lim^^oo '''j = 00. Due to p.24p . 
the right-hand side of the last equation in p.ip is bounded. Therefore, there exists a subsequence 
of {tj}, that we label in the same manner for simplicity, such that: 

inf |*(t)| ^ 00 . 

This means that infT-£[^^_i ,-^.+1] -^^ — > 00. Therefore sup^gr^._]^^.^]^i $(t) — > 00. Indeed, if 
$(tj ± 1) are bounded, and since we can determine uniquely $(t) in the interval re (tj — 1, tj + 1) 
using $(rj ± 1) and ^^{t), t S [tj ~ 1, tj + 1] the claim follows. This contradicts p.24[) . whence 
limsupT._>g2 l^(''')l < 00. Then, (|3.25p implies lim^_i.oo W^(''') = 00. Using the first equation in 
()3.ip it then follows that ^(t) is unbounded, and this contradicts again p.24p . whence (|3.15p 
follows. 

It only remains to prove p.l6p . We first notice that 

lim inf $(r) < 2 . (3.26) 

Indeed; otherwise it follows from the last two equations of p.ip that ^^ < ~^! this implies that 
$(r) vanishes for a finite value of r. This contradicts Lemma [5J 

We now claim that limsup^_j.g^ $ > 1. We argue by contradiction. Suppose that 

lim sup $ < 1 , (3.27) 

then, using the last equation in p.ip . we have ^ > (5 > for r large enough. This means that 
lim^_s.oo ^(t) = c». Using again p.ip . we obtain limT-_j.oo W{t) — cxj and lim^_yQo ^{t) = 00, this 
contradicts (|3.27l) and therefore implies limsup^^^^ $ > 1. 

We now claim that lim sup^_j.o2 $ > lim inf t- _j.oo ^ (''')• Suppose to the contrary that lim^_j.oo ^(t) 
^00 > 1 is defined. If £00 > 1, we obtain a contradiction because p.ip . in that case, implies that 
$(r) vanishes at some finite r, but this contradicts Lemma [2] And the case £00 = 1 is against the 
assumption of the lemma. 

We have then limsup.^_^3^ ^{'t) > lim inf ^^yoo ^{t)- This implies the existence of sequences of 
minima and maxima of $, namely, {t+}, {t^} and {e,i}, such that lim„_j.oo t+ ^ lim„_j.oo Tn — ^^j 
£„ > for all n, and such that M^(t+) = W{t^) ~ with max^^r +_^ ^+^^ 1 $(r) = $(t+), 

min^gT +_^^ ,^+_^^^i ^(r) = $(t^) and lim„_^oo *(t+) = limsup.^^^ $(t), lim„_^oo *(t^) = 

lim inf T-^oo '!'('''). Suppose that lim inf t-_j.oo *!'(''') > 0. Due to p.26p and the definition of E 
in p.l7p it then follows that 

1 ,, . 1 



Eir-) = ———-r + HT-)< 



2(<i>(T-))2 '"'- [liminf,^oo$(r)]' 

for n sufficiently large. This contradicts that E{oo) = cxj and the second formula in (j3.16p 
follows. Suppose now that limsup^_j.g^ $(r) < 00. Using again the definition of E and that 
lim sup^_j.g^ $ > 1 we obtain 

^("■» ) = o^^/ + NN2 + '^(^n) < 1 + 2 hm sup <J>(r) 
which also contradicts that E{oo) = cx). ■ 
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3.2 Oscillations: The formal description 

We now formally describe in some detail the oscillatory solutions of (|1.3p . Such solutions have 
increasing amplitude of oscillation as r increases (see Lemma [5]) . When this amplitude becomes 
rather large as r — >■ oo, one expects the behaviour of the solution to be dominated by (|1.5p . This 
equation can be solved explicitly: 

$(r) = -TV6 + ciT2+c2r. (3.28) 

Thus the solution in such regimes resembles a third order polynomial near a local maximum. For 
larger values of t the magnitude of the solution becomes small and its behaviour is dominated 
instead by (|1.6|) . but, as it is proved in Lemma [21 the solution does not vanish for finite values 
of r and it then increases back again to a larger amplitude. How this happens is described by 
(|1.6|) . Thus the function $(t) can be described by means of alternating regimes, where either 
<I>(r) becomes very large or $(t) is close to zero. As it will be seen later, the matching between 
such regimes will require that, when $ becomes small for increasing r, the polynomial of the form 
p.28|) that is asymptotically close to the solution, has a simple zero. On the other hand, when the 
solution leaves a local minimum (sufficiently close to for r large enough) the polynomial (|3.28p 
that it approaches asymptotically has a double zero. These, linear and quadratic, behaviours near 
$ ^ result from the analysis of (|1.6p . this is done in Appendix 2, see Theorem [T5l 

Since $(t) is oscillatory (see ()3.16|) ). there exists an increasing sequence {r+} of local maxima, 
thus having ^{t^) = 0. For simplicity of notation, for each n, we define 

Then Lemma [5] (|3.16p implies limji_j.oo ^n ~ oo. As indicated above, $(t) must have a double 
zero, at least to the leading order as L„ ^- oo. Using the fact that $ solves p.5p to the leading 
order, it can be seen by means of an algebraic computation, that the only possibility is to have 
the following asymptotics if L„ — )■ cx): 

, t2 

^{t+) = L„ , ^<J>(r+) = , ^*(r+) = -aLl , (3.29) 

since $(t) must have a double zero at some t < t^ , to leading order as Ln — > oo, a will be chosen 
later in order to fulfil this condition. We introduce the scaling, suggested by (|3.29p . 

$(t) = Ln^niVn) , Vn = ^^^ (3.30) 



where ^niVn) satisfies 



that, to leading order, becomes 



4*- + '^l3(i^' <^-^'' 



dVn 



3 $„ + 1 = , (3.32) 



with initial conditions 



$„(0) = 1, ^(0)=0, ^(0) = -a. (3.33) 



Solving (|3.32p - (|3.33p gives the approximation $„(?/„) — Poo('7n), where 

P^{v)^~lv'~^v' + l- (3.34) 

Imposing now that Pqo has a double zero at a negative rjo requires: 

«'n(?7o) = - g»7o -^^0 + 1 = 0, Kim) = -^Vo - a7?o = , 
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770 = -2a, a=(^) • (3.35) 



On the other hand, the polynomial Pooij]) vanishes at a positive value of rj, namely at r]a 



:ir- 



We can then approximate the function $(t) in the intervals where t G ( t+ + rjoLn, t+ + rjaL. 

by means of (|3.30p and p.34p with p.SSp . The approximations of § by (|3.34p near ri~ and near 
rjQ arc then as follows 

d - /3\ ^ 

*,i(?7n) K{T]n-r]a) as r]n^r]- with ^ = "j— *"(^") = ( 9 ) ^^ ^^'^^^ 



and 



a , ,9 I. (P - , , /3 



$«(??«) ^ -^(j?™ -?7o)^ as r]n^T]^, since -^^^"(r/o) -^ ,„ 
2 d774 V2 



In order to describe the function $(t) for the values of r where $ becomes small we introduce 
an inner layer variable near rja for every n. Using p.36p we can infer that this inner layer is 
characterised by 

Vn = Va + YjCn ^{Vn) = Jj'finiCn) (3.37) 

with (fiiCn) satisfying, to leading order if n is large enough, the equation 



with the matching condition 



(3.38) 



^n{Cn)--KCn+^-^j0^ + ^ + o{l) aS Cn ^ -^ , (3.39) 

where the logarithmic correction and the next order has been obtained, computing the next order 
in the asymptotics of the solution of p.3ip and p.33p . More precisely, this correction is given by 
using the leading behaviour p.36p . namely 

rVn i-si r^'^ I 1 3 

/ / / 7^— 73^*3rfS2rfsi = --— 3l0g(|77„-77a|) + -— 3 +0(1) aS rjn^Va 

Jo JQ Jo ($„(S3)) ^^ ^^ 

where $n(s3) is as in p.34p . 

Using Theorem [T5] in Appendix 2, it follows that there exists a unique solution of (|3.38p . (|3.39p 

given by (pniCn) = </? Kn - '"2^4" - 5^ j . Therefore: 

<^„(Cr0^rC'. as Cn^oo, (3.40) 

where F > would in fact be fixed by the matching region. Wc observe that the matching 
condition (|3.33p is valid for |^„| ^ log(i„), |Cn| ^ ifi, Cn < 0. On the other hand, the matching 
condition p.40p is valid for |C„| ^ log(L„), (n > 0. In order to determine the maximal region 
of validity for this asymptotics we need to find the size of the region where the second term on 
the right-hand side of p.3ip becomes relevant. A standard dominated balance argument indicates 
that this happens for r]n — r]a ~ 0{L'^J. We then use the following change of variables: 

^n{Vn)^(^) Ll^n+liVn+l) , Vn+1 + 2a = ^^^L^ (3.41) 
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where a is as in (|3.35p . The particular rescahng has been chosen in order to have <i>„+i(0) = 1 to 
the leading order. Notice that, using (|3.3ip we obtain: 



(P - _ / a \9 1 1 



,3 <^n+i + l=i^ 7^7— -3 (3-42) 



Using (|3.37p and (|3.40p we obtain the matching condition for $n+i(^/n+i): 

^n+il^yn+i) ^ 2 ('^"+1 + ^") ^^ '^"+1 ~^ ~'^"' (3.43) 

and the solution of (|3.42p . (|3.43|) to the leading order gives the approximation ^n+iiVn+i) = 
Poo{Vn+i)- Notice that we obtain $„_|.i(0) = 1 to the leading order, as expected. Note also that 
(lOIl . (|3^ imply: 

$(r) = f ^y L,io$„+i(rM+i) , r7„+i = -2a ' ° ^^ " ^^^ " '''^^" 



2r 

Equation p.42p is equivalent to p.3ip and therefore we have 



10 

T 3 
J-in 



^n+i=(f)'c (3.44) 

On the other hand, since ?yn+i ~ at r = t^j^^ we obtain: 

T++1 = T+ + r^ail + 4rL? (3.45) 

The sequences (|3.44p and p.45p describe how the increase in the amplitude of $ takes place. 
Notice that: 

L^=^^^^^^^{l + o[l)), r+=2af(2r)iexp(§(10)")(l + o(l)) as n ^ oo 

(3.46) 
This formula can be easily obtained from p.44p using the fact that Zn = log(L„) solves the 
linear recursive equation 

Z,,+, = 10Z,, + 3\ognpj . 

Moreover, suppose that we denote by Af„ the minimum values of $(t) in the interval {t^ , t^-^^^) , 
and that these minima are reached at the points t^ G (tJ",t^+i)- Then, using p.30p and p.37p . 
as well as the arguments yielding (|3.45l) we obtain: 



mina ip{C) , 



M^ = t^^^^^^P^il + oil)). 



2a^(2r)f ■ ^'^'^^ 



exp ( — (10)" ) (1 + o(1)) as n -^ oo 



i2T)i 

(3.47) 
where f is as in Theorem 1151 

3.3 Oscillations: The rigorous construction 

In this Section we prove that every solution of p. 31) such that $(t) docs not converge to 1 as 
T — > oo oscillates with increasing amplitude as r — > oo in the form described above by the formal 
asymptotics. Namely, we prove the following theorem: 
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Theorem 6 Suppose that $(t) is a solution of \1.3\) that is not included in the stable manifold 
of the point Pg. There exists two increasing sequences {t^}, {t^} satisfying t^ < t^ < t^^-^ and 



lini„_>oo ■''n ~ °°' such that: 



d^ 

— (r) > m (t„ , T+) , $, (t) < m 



Vn ' '''n+l) 



Moreover, for every n, let i„ and A/„ be 

L„:=$(r+), A/r„ := <J>(t-) 
Then, the asymptotics \3.4^()^ and \3.41^ hold, and also: 

^^T^^-U) m^ + o{l)) as n^oo 



(3.48) 



(3.49) 



Remark 7 Notice that Theorem [^ implies that every solution of lil.3\) that does not approach 
(f> = 1 as T —> oo oscillates for large values of t in the precise manner indicated in the Theorem. 

The proof of Theorem [5] will be decomposed in a series of Lemmas. We first observe that the 
existence of the sequences {t,^} and {r,^} is a direct consequence of Lemma [5l Also, if L„ and 
M„ are as in (j3.48p for all n, then Lemma [5] implies that lim„_j.oo Ln = oo and lim„_j.oo Mn = 0. 
We define the sequence {a„} C M by means of 

-TT^iTn) = -anLn , whcrc a„ > 



and we shall then prove that an ^i' a, with a as in (|3.35p . as n — > oo. This then implies p.49p . 

We define a sequence of functions ^n{Vn) as in (|3.30p . but where $„ solves (I3.3ip with the 
initial conditions: 

d$r, ,„, . d'^<^r 



$„(0) = 1 



-(o)-o, 



-(0) 



(3.50) 



drjn^"' "' dvil 
Observe that the only difference with p.33p is that the last condition depends on n. We define 

1„3 _ a^ 
6 ' 2 



the approximating polynomials and their roots accordingly; letting Pa^{ri) = — ^?7^ — ^77^ + 1, we 



define r]a„ as the solution of 

We observe that when a„ is very large then rja^ ^ 1/ y/a^, and otherwise ria„ is of order one. For 
that reason and in order to avoid studying several possible regimes for a„, the scalings with L„ 
would include a factor depending on On, as we shall see below. 

We first derive some approximation formulae for the functions ^nivn) in the region where they 
are not too small. 

Lemma 8 Let $ a solution of il.3]} satisfying the assumption of Theorem\^ Let $,i(7]„) be as in 

) independent of n and of a„ 



I13.30\) . Then, there exists a N > independent of n and of an such that the estimates 

1 



2?7a„ 



^n('7„) - 1 - 



d<^n{Vn) 



anVl 



6 



< 



drin 

d'^^niVn) 



anVn + -^ 



dVn 



< 



< 



4(^aJ=^ 

LI 



log(l-^ 



1 



Ll iVa„ - ?M) 
LI {Va„ --HnY 



hold as long as: 



- L3 



l°g (1 - ^) 



(3.51) 
(3.52) 

(3.53) 

(3.54) 

(3.55) 
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Proof. Direct integration of (|3.3ip with p.SOp gives: 

I rVn /-Sl /•S2 



drjn 

dVn 



2 3 

^nVn Vn 



1 



2 6 ' LIJ, J, J, ($„(S3))3 

,2 j^ rVn rsi ^ 



ds3 ds2 dsi , 



^nVri 



Vn 



^n ^?n T J- o 



^« 7o Jo (^n(s2))' 



-ds2 ds 



1 J 



(3.56) 



i;i Jo (*«(si))^ 



We now claim that p.5ip - p.54p hold if p.55p is satisfied for some iV > independent of L„ and of 
a„. Indeed, this is proved by means of a continuation argument. The inequality (j3.5ip is satisfied 
for rjn = 0. On the other hand, as long as this inequality is satisfied we have: 



^niVn) - 1 



an?7„ 



6 



< 



LI 



■rin fSl /•S2 



Jo Jo 



1 



iVa„ - Ss)^ 



ds3 ds2 dsi 



and hence p.52p also follows. Due to the convexity of the polynomial Pa„{v) we have Pa„{rin) > 
■ncL„ 



It then follows that: 



^n{Vn) > 



iVa„ - Vn) 4(77a„)^ 



Va„ 



Ll 



log 1 



Vn_ 
Va^ 



which implies p.5ip . Moreover, using p.56p . we obtain that p.52p . ()3.53p and p.54p hold as long 
as (|335l) holds. ■ 

It is now convenient to reformulate this result using a new set of variables, which is better 
suited for the study of the boundary layer where the term -p- , . ,3 in p.3ip becomes relevant. 
Namely we take: 



^n{rin) = -p^^n\(,n) , Vn - Va„ = r^ C» 



and (|3.3ip becomes 



d^ifr, 



Ll 



Ll 



(3.57) 



(3.58) 



da Ll ivnr 

The next step is to change into the variables which transform (|1.6p into a planar system of 
ODEs. We use the transformation (JA2.5P of Appendix 2, that now reads 



dfn 
dCn 



fn Un , 



da 



Vn Un : 



C« 



(¥'n(s))3 



where 



RM,n = M 



loe 



Ln 



for some Af > 1 , 



(3.59) 



(3.60) 



(in this way, C„ very negative is in the matching region, see p.55p ). This transforms p.58p into 
the system 



d^pn dun 

— = Unfn , "1 = V„ 

dZfi fJ'Zfi 



1 



dVn 
dZn 



= 1 + TiUn , Vn 



Ll 



{^nf , 



(3.61) 



that can be considered as a perturbation of the system (JA2.6P for n large enough. We will use 
u = Un SiS independent variable and one should understand in the following, with some abuse of 
notation, that (p„ = ipn{u), v„ = v,i(u), Zn = Zn{u), then 



dcpn 


Uifin 


dVn 

du 


Vn + |u^ 


-(V'n)^ 


dZn 


du 


Vn + 1^2 ' 


du 



1 



Vn + ^U'^ 

The following lemma is a translation of Lemma [8] into this system in the matching region. 



(3.62) 
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Lemma 9 (i) Suppose that (u„ + iu^j) is large. Then we can define (y9„(u), u„(w), Zn{u) by 
means of \3. 59j) with u ~ Un- 

(ii) For any £ > arbitrarily small, there exist Mq ~ M(){e) and uq = no{e,M) such that for 
M > Mq, n > uq we have: 



\(p„{u*) - R'nRM.nl < sRaLu , \Vniu*)\ < £ , Zr, 







(3.63) 



for some u* £ 



-i^|(l + £)(i?A/,„) 3 , -i^|(l - £)(i?A/,„)5 



with 



Kn= [ anT]a„ + — 



Proof. The possibility of defining the functions fn{u), Vn{u), z„(u) is just a consequence of the 
second equation in p.61|) and the Imphcit Function Theorem. Using Lemma [5] we obtain: 



Wn 



and 



KnCnl <e|Cn| , 
d'^ipn 



dC,n 



"■nVL 



da \ LI 
as long as n is sufficiently large and 



Kr. 



Ll 



< C{a„vlJ 



if 



Ln 



\Cn\ 



4 



ICnl 



C^ > ICnl > M 
Van 



log 



Ln 



where C > is independent of n (the inequality on the left-hand side is used in the first and 
second inequalities). Moreover, using that the sequence a„?7^ is bounded, as well as the fact that 
Va„/Ln — > as n — !• oo, wc obtain that at Cn = ~Lti\i.n 



\ipn - K„RM,n\ < sRm,, 



dQn 



Kn 



< £ and 



d'^(Pn 



dQ 



< 



Ll 



7?2 



hold, where e can be made arbitrarily small if we choose A'l large enough and n > no{M) sufficiently 
large too. Thus, the functions {^n,UmVn) satisfy 

|(/3„ - KnRALnl < sRhLn , "n + K^ {RhLn) " < £(i?A/,n) " , \Vn\ < £ at Cn = ~RM,n 

(3.64) 
where e can be made arbitrarily small if M and n are large enough. Defining u* as the value of 
Un at which C„ = —Rmji we obtain p.63p . ■ 

We now need to prove a stability lemma for the system p.62p . The behaviour of ((/?„, u„, w„) is 
expected to be similar to that of the separatrix parametrised by {(it, v) : v = v{u)} of Lemma [TBI 
of Appendix 2 and that gives the linear and parabolic behaviours in terms of ip, see Theorem 1151 
Before we go into such analysis, we prove some results for the following auxiliary functions. 

Lemma 10 Let us define functions G and H by means of 



G{u) = exp 



Jo V{s) , 3 

Then, the following asymptotics hold 



-ds 



H{u) = exp 




G{u)=K+{uy^\l + 0\^ 



G{u)=K_\uf il + o' ^ 



as M — > oo and 



as u — > — oo . 



(3.65) 

(3.66) 

(3.67) 
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with Eq > and where K+, A- > 0. We have also that G{u) > for all u Cz M. and we can write 
the following global estimate: 

■^ (l + \uf^"'') < Giu) < C, (l + luf^'^A , /3(w) = - ifu>0, /3(w)==3 tf u < 

(3.68) 
where C* > 1. Moreover, 

H{u) — K+{u)'^' I 1 + O I \\ as M — > oo and 

H{u) = K^\uf(l + o(-^)\ as M^-oo, 

with Eq > and where K^, K^ > 0. And, as above, there exists a C* > 1 such that 

^(l + l^r^"^) <if(u)<a(l + |ur^"^) , J{u) = ^ if u>Q, 7(w) = 5 if u<0. 

(3.69) 

Proof. This result is just a consequence of the asymptotics of the function v — v{u) obtained in 
the Appendix 2, see Lemma [T51 (|A2.7p . ■ 

Lemma 11 Suppose that (fn{u) and Vn{u) are defined as in Lemma \^ Then, for any e > 
there exists Mq = Mo(s) and no = no{s,M) such that, for M > Mq and n > no the following 
inequalities hold: 

\{vn-v)iu)\< , ^^^ ,x g(»), (3.70) 

i + iu*r 

with Ci > large enough and independent of n, e and M , and 

ip„{u) G{u) 



ifniu*) G{u*) 



< ^^^#^ (3.71) 



[l + \u*\')G{u*) 



1 f r-™ 

for u G [u* ,u**] with u** ~ 5~5^ '^^'^ u* *■' '^■5 ^'^ Lemmai^ and C > is independent of e, n, M 
and Ci . 

Proof. Using ([3:62|) and (|A2.10|) we obtain 

ip„{u) = ip„{u*)cxp\ / — — ^rds] (3.72) 

\Jw Vn{s) + % J 

and 

d{v^ -«)_(-! + ¥) («" - ^) (r?,„ )9 (^„)3 



Then, p.7ip follows from p.72p and Lemma [TUl (|3.69p (using Taylor) as long as p.70p is satisfied. 
Let us now prove (j3.70p . First observe that we can write 



du 
with 



^^--^^)^ti^^%^./(.) (3.73) 



f^^^^t'^'^)^""'"U 1 1 ^ iVaS {^n) 
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Thus, integrating p.73|) we obtain 



-^ ^ 9 



IK -C)(w)| < 2eexp I / ^ ^ds | + / exp | / ^ ^ds | f{a)da. (3.74) 

H^) + t) 

Also, using p. 651) and (|3.69p . we have 




5s-^ 



exp / ^ ^-^ds < , \ ■ (3.75) 

Notice that p.70p holds for u ^ u* because Lemmal^ p.63p implies that \vn —v\ < e aX u = u* 
for n sufficiently large. Therefore, we can extend it to some interval contained in u > u* by means 
of a continuation argument. In particular, if p.70p is satisfied then also 

, 1 / u" 

K--v\<^\v + - 



is satisfied, and hence 



\f{u)\< 



9 fc^.^3 2(1 + ^ (z;„ 



LI {v + f 



Then, using (|3.7ip and as long as p.70p holds we can further estimate f{u) as 

\J{U)\ S T7. ^^/ ,*NN,C( r- . „2\ + ^f^l^ 



L6 (G(^*))3 (,, + i^) ^ (z; + f)' (l + kl^°)' 

Using this and p.75p in p.74p we obtain 

2£aH(^) 2(r;,j9 (y,„(^*))3 /•« iJ(^,) (G(a))3 ^ 

(Wn - w (u) < -— — -^ + — — / TTTTT^T"^^^^^" 

(1 + |u*|5) Lb (G(u*))-^ Ju' H{a) (v{a) + ^) 



' A. i/(a)(^(^) + ^)3(l + kP°) ^ ' 

Now, using (|3:68| and ([3:69| in ([3J6| gives 

.,, „^ 2C*£g(») , ^ (7?aJ^(y„(u*))^ ^, ,,^ , , „3, , 

IK - -)HI < (Tip^-rpy + ^-^TgRF ^"^ ^ + '" ' ^ ^ 

^ LI (G(u*))3^("H^ + ("M+^^i" (l + KD ^^-^^^ 

where (it)+ = u \i u > i) and (u)+ = otherwise. Notice that for u < 0, the estimate p.77p 
implies p.70p if we assume that Ci > 3C* and n is large and e is small. In this case the second 
and third terms in the right-hand side can be absorbed in the first by noticing that ^^ behaves 
like exp(— i3(u*)'^) for a positive constant B and that (^„(m*) behaves like {u*Y (see Lemma[9|, 
so, in particular 

Ll {G{u*)Y ^ l + |u*|8 ^ l + |u*|5 

for n large enough, and where we also use p.68p . Therefore, a continuation argument implies that 
p.70[) for u < 0. In order to derive the range of values of u > we only need to see the range of 
values for which the third term in p.77p can be absorbed into the first one, that is 

^^(u)J < C- — ^ with C> . 

L6 ^ '+ 1 + |m*|5 
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Thus taking for example, u < —^ for n large enough the estimates holds (absorbing a factor Ln 

in order to control e factor). ■ 

In terms of the variable u and the transformation p.59p , the balance that gives the maximum 
size of the next matching region corresponds to the one for which the second and third term in the 
equation for ^^ in p.62p are of the same order (that is |uu„ ~ 2^ (y'n)'^), thus u* < u <^ -^ 

as n — > oo. On the other hand u** = ^ , ",5 < 7 — ^W for n large enough. 

As a next step we obtain approximations for the function $,1(77,1) and the derivatives — nilhil 



dri„ ' 

— ^" i" to the most right value of rjn in the boundary layer region where the function $„ is small. 
More precisely, this most right value of r]n will be defined by means oi u** . 

- / \ I 
Lemma 12 Let ^^,^, ~ -^ ( -tt^ ) (m**)5 . For any S > 0, there exists no = 'no{d) such that for 

n > no{6), there exist 77** £ iia„ + ta (1 — <^), ?7a„ + A (1 + S) such that: 









< 



< 



dVn 



2(^a„)^ 



< 





2 









(ryaj^ 



(3.78) 



where A„ 



m' 



K^ 



Remark 13 The coefficient T in Theorem \15\ can be approximated from ^3. 7<$| j for n large enough 
as T ^ A„/4. In particular this is consistent with relation T ex K^ (which follows by a scaling 
argument in \A2.1\) and the behaviours iA2.S\) and IIA2.3\) . 



Proof. Notice that we can assume that z„ in p.59p . and then also Cn are functions of u for 
u* < u < u** . Combining p.59p . (|3.62p and using also (^„ — ^Rm.u aX u — u* we obtain: 



Cn + RhLr. 



w w„(s) + ^ 



-ds . 



Clearly, from p.63p and (|3.67p . '^/".^^ is uniformly bounded by some constant C > 0. Then, using 
Lemma ITT] (estimates (|3J0| and (pTTTj) ) and (|A2.7j) . 

'Vn{u*)\-^ r {Gis))i 



Cn + RM,r, 

CC\e 



< 



Giu*) 



■ds + 



v{s) + f 
CCie 



-ds 



iGis))lHis) 



ds for u < 



il+\u*\^)J^, 1 + ^2 • ■ (1+K|5)7„. (1 + ^2)2 

Using also the global estimates (|3.68p and p.69p of Lemma [TOl we further get that 



(3.79) 

(3.80) 



Cn + Rm.t. 



Mu*)Y r (Gis))! ^^ fM^*)Y r (Gis))i 



G{u*) J J^,v{s) + ^ \G{u*)J Jo v{s) + ^^ 



ds 



<CCie{l + \u*\) + 



CCie 



l + {u)l] for u<u** 



(1 + KP) 
where we have split the integral on the left hand side of p.80p . 



(3.81) 
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Now, using dMl), (EMI), (P:B7|) and ((Mj)) . we obtain: 



«„„.-,^3Mvri2(fe.» 



G{u*) J J,,.vis) + f 



< SRM,n 



(3.82) 



where (5 > can be made arbitrarily small if e is small and n is large. On the other hand, using 
(|3J4l) . ([3:661 and (|A2J| . we obtain: 






G(w*); Vo ^;(s) + f K\K. 



< -u5 u > 0, 
- 2 



(3.83) 



where 5 can be made small choosing e small and n large. Combining now p.8ip . p.82p and p.83p 
we obtain, for small e and large n: 



<5(RM,n + {u)-']+CCieil + \u*\) for u<\u**\ 



^""i|(l^)'^(")* 



In particular, recalling that u** ^ cxpB{u*)^ for some positive B, this gives 

2 /K^ 



Cn 



K^ \K 



{u**Y 



< 6{u**)^ 



(3.84) 



with S small if s is small and n large. 

On the other hand, using p.7ip . p.66p . and similarly using ()A2.7p . p.69p and p.70p . we obtain 



Vn (U ) - 



K+ iu**)i 



K- Kl 



< (5(u**)5 and 



Vniu**) 



**^ iu**r 



<6{u**y' 



(3.85) 



with 5 small for s small and n large. Translating these behaviours by means of p.59p . we finally 
obtain 



dC,n 



K_ 



A'„(m**)5 



<5{u**y^, 



dVn 1 fK_ 



da 2 \K+ 



Kl 



< 6 at u = u* 



These estimates, as well as p.85p . can be written in terms of the variable C,n by using p.84p . 
namely, 



at some C„ € 



,iK-Ya 



K+J 4 

2(] 



< sC 



d(Pn r^^5 f^-\'' Cri 



dCn " V^' 



< SCn 



dVn 1 fK- 



dQ 2 yK 



K' 



<S 



Observe that i^„ can be bounded 



(^(^)^.-)i,HM£)(|±)^.-)^ 

above and below by constants or order one, independently on the choice of a„. The result (j3.78p 
follows now from p.57p . ■ 

We now approximate $«(??«) for rjn > rj**. To this end we define the following polynomials 
that approximate, as we shall see, the functions $„ for n large enough in the appropriate intervals. 
Thus we let 



Qn{Y) 



4 " 2 ^4 6 



Bn = ^iVa„f, 



where An is as in Lemma [T2] The definition of ^** (cf. Lemma [T2|) as well as (|3.60p and (|3.63p 
imply that lim,n.oo Bn ~ 0. The next result now follows: 
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Lemma 14 For S > 0, there exists no{S) such that for n > no{S) the function ^niVn) is increasing 



in an interval rjn G (??,**, vT^Vn), where rjn G 
Moreover, defining in+i = ^{Tn+i) '^^ obtain: 



^)'«)*i^(i-*).(fe)' 



^2^ (')a„) 



.(l + <5) 



in4 









(3.86) 



where y*^„ = i I ^„ + \/ A^^ + 4A„ i3„ j , and afeo; 



'n+1 'n '/ti ^n ^*,n 



{ria^f 



< SY., 



Lrf 



(3.87) 



Proof. First, we define J„ — 
Integrating (|3.3ip we obtain 



i^?::7T 



:, for simplicity of notation, and observe that lini„_^oo Jn = oo. 



'^niVn) =PiVn) + 



Ll 7^.. 7^-. 7^.. (^n(s3)) 



(is3 ds2 dsi 



(3.88) 



with. 



d7]n 2 dr/f^ 



6 



As a direct consequence of Lemma [T^] one has that 



|^(^n)-^'Q(>^)| 



d7]r, 



(PiVn) - JlQ{Y)) 



dv?i 



(PiVn) - JlQ{Y)) 



< 5JlQ{Y) (3.89) 



with Tj** < rjn, Y < (1 + 5)y*.„, if n is sufficiently large and where r/„ — r/** = Jnl^. 

We can now estimate the effect of the last term in p.SSp . To this end we use the type of 
continuation argument that we had used repeatedly. Now our goal is to show that: 



$n(??n) > 



J?.QiY) 

8 



(3.90) 



This inequality is satisfied for rjn — rj** ■ On the other hand, as long as this inequality holds, we 
have, using the change of variables rjn ~ rj^* ~ JnY , that 



Ll 7,,-* 7^-* 7^;- (*n(s3))^ 



1 29 

-ds3 ds2 dsi < 



Y rYi rY2 



LlJUo 7o 7o (QiYsW 



rdYs dY2 dYi 



If < Fa < >;.n(l + S) then QiYs) > C (^ B^ + \ B^Y-^ + ^Y^). Then, using the change of 
variable Y = B^Z we obtain: 



1 



1 



< 



C 



^n'^rfin JO 



y/B„ fZi 



Ll 7,,** Jv'„' Jvr (*"(s3)) 



dZ^ dZ2 dZi < 



,T + i-^3T- T^3 



Ay2\^ 



'-ZI 



-dsz ds2 dsi 

r3 76 d5 ' 



for n large enough. Here the first integral in the integrand can be estimated by a constant using 
that Bn — >■ as n ^ oo. Analogous estimates can be proved for the first two derivatives of this 
integral with a similar argument. 
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Wc observe that ^^^^g " , ^ J^ as n -^ oo, since this is equivalent to C ^ L^J'^B^ ~ (C**)^ 
as n ^ oo and this is satisfied due to the definition of ^»*. Therefore p.90p holds for < ^ < 
(1 + (5)i;,„. It then follows from ([XSS)) . ((XMI that: 



\^n{71n) ^ JIQ{Y)\ 



-^i^niVn) - 4Q{Y)) 



^(MVn) - J^Q{Y)) <dJ^Q{Y), (3.91) 



with r/** < rjn and F < (1 + (5)y*,„. 

It then follows the existence of 77„ > as in the statement of the Lemma, and such that 
rf|:(^nte* +fjn)) = 0, and that ^($„(7y„)) > for rj^ G «*,<* + fjn). Moreover, dSH]) 
implies (|3^ and (jXSOl) yields ([3:87)) ! ■ 

Proof of Theorem [6j The asymptotics p.46p is just a consequence of p.86p that can be solved 
recursively arguing as in the formal derivation of (|3.46p . On the other hand (|3.47p is a consequence 
of p.30p . p.57p and Lemma fTTl Finally p.49p is a consequence of p.9ip . Indeed, notice that this 
last formula implies 

ar.+,=-{l + 0{S)) ^'^*l , n,. = ^(A„ + v/^2+4^^5j. (3.92) 

(Q(i;))-^ ^ 

Since Bn — > 0, we have (y*,n — A„) — > as n — )- oo. Therefore < eo < On+i < — for ti large. 
Then 

g"(n,„) ^3^ ' 

"" " ' oo 



(Q(n,„))^ V2. 

whence (|3.92p implies, since S can be assumed to be arbitrarily small for large n, that an+i — >■ a, 
with a as in (|3.35l) . and the result follows. ■ 
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Appendix 1. The operator (VF + VF^)N in curvilinear coor- 
dinates 

As above, vectors components in Cartesian coordinates are denoted by sub-indexes a; and y respec- 
tively, and when expressed in curvilinear coordinates the indexes 1 (tangential to the substrate) 
and 2 (normal to the substrate) are used. For a general vector F, thus with F = F^e^ + FyGy = 
Fit-|-i^2n, we compute (VF + VF-^)N. This is done, for example, using the expression in curvilin- 
ear coordinates for linear differential operators given in e.g. [1], and one obtains, writing separately 
the normal and tangential components: 

((VF + VF^)N) •N = 

2g^, f^i-a~;f. (^JL.kF^+2il-kh)'^^'^^^'''^^'^'^^ 



([i-khY + {^f)li-kh)\^^ J {i-khr 



.ds, 



2f^ii-kk)uHm -ji-khn. (§+.F,+(i-fc.)^ 



{i-khr + {^y){i-kh) 
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and 



((VF + VF^)N)-T 



>f (l-fcfe)^, + 2((l-fcfe)^-(f )^)ti fdF, , ^^.^^ ^ ^^ ,^^9Fi 



((l-fc/i)2 + (f)^)(l-fc/.) \ds 



Appendix 2. Analysis of the equation (11.61) 



In this section we concentrate on the study of equation p.6p . In order to clarify them matching 
conditions needed in the inner regions where $ ~ that are part of the oscillatory solutions to 
(jl.3p . The inner variables used in this inner regions are Cn, and (pn and so we rewrite (|1.6p in this 
variables but dropping the subindex n. The main result of this appendix is: 



(A2.1) 



Theorem 15 There exists a unique solution of 

d^ip _ 1 

with the matching condition: 

ifiC)'^ -KC + o{l) as C^-oo. (A2.2) 

Moreover, the asymptotics of (^(C) for large C is given by: 

ip(Q ^ rC^ as C ^ oo for some T > . (A2.3) 

Finally, the exists a unique solution of \A2.9\) with matching condition 

vie) ^ ^KC + o(l) o,^ C, ^ GO for some K > Q . 
It also satisfies that there exists a finite Q* such that 

^(O^O as C^(C)+. (A2.4) 

All other solutions satisfy ^A2.!^) for increasing C,, and, for decreasing C, either IA2.4^ or 
(/9(C) ^TC, as ( —^ — oo for some T > holds. 

Wc remark that in this paper we only need the existence of the solution satisfying (|A2.2[) and 
(|A2.3|) , the rest of the behaviours are stated for completeness and because they will be used in [9] . 

In order to prove this theorem, we reduce equation (|A2.ip to a second order autonomous system 
of ODEs, as follows 

— = ip=>u,—^ip=>v,dC^^p^dz, (A2.5) 

giving the system 

d(f du 1 2 dv . ^ , ,c^ r.\ 

_=u^, _=, + _, , ^ = l+3«u. (A2.6) 

We now prove several lemmas that give the behaviour of the orbits of this system. 

The last two equations in (|A2.6p can be studied independently by means of a phase-plane analy- 
sis. The isoclines are Fi = {(ujii) : u-t-iu^ = 0}, that has dw/dz = 0, and r2 = {(u,?;) : l + ^vu = 
0}, that has dv/dz = 0. The only critical point is p^ ~ {u^,Ve) — ((9/5)3 , — (l/3)(9/5)3 ), and lin- 
earisation gives two complex eigenvalues with positive real part, namely A± = g (|) ^ (7 ± vTli). 
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Figure 1: Phase portrait associated to (jA2.6[) showing the direction field. The sohd hnes represent 
the isochnes and the dashed ones the separatrices. 



Wc distinguish five regions in the phase-plane, the ones separated by the isoclines. These 
arc depicted in Figure [1] where direction field is also shown in each of these regions. Standard 
arguments imply that any orbit on the phase-plane eventually crosses the isoclines into the region 
i?4 forwardly in z. If, however, an orbit has (m, v) e i?4 at some value of z, it is possible to discern 
from which of the regions is coming from for smaller values of z by identifying the separatrices of 
the system. We have the following result. 

Lemma 16 (Separatrices) (i) There exists a unique orbit v ~ v{u) in the phase-plane as- 
sociated to system HA 2.6]) that is contained in R4 for all u G M. Moreover, v{u) has the 
following asymptotic behaviour 



v{u) 



O u 



as u 



00 . 



v(u) 



2u 



(1 + 0(1)) as u 



-00 



(A2.7) 



(ii) There exists a unique orbit v = i){u) in the phase-plane associated to system liA2.6\) that has 
the following asymptotic behaviour 



v(u) = — -— (1 + o (1)) as u 
2u 



+00 



and 



|(u,w)-(ue,We)| <Ce^^(^+)^ as z^-00. 
Remark 17 The separatrix v ~ v{u) is uniquely determined by the problem 



dv 
du 



v{u) ~ — — as u 



V + iu^ ' - \ ' 2u 



(A2.8) 
(A2.9) 

(A2.10) 



Proof. In order to capture the behaviour (jA2.8|) . wc perform the following transformation. 



w ~ uv , — - = V—- 
dr dz 



(A2.11) 



If the behaviour is as expected then the transformation is locally valid as long as w 7^ 0, thus, one 
can integrate the second equation in (IA2.11[) . Using the equation for u in (JA2.6P we obtain the 
system 

--^w + 2w^ -{-v^ , —=v[l + -w], (A2.12) 

dr dr \ 3 J 
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that has critical points (— ^,0), (0,0) and (— |,— (^) ^ )• The point (O, — ^) is a saddle-point. 

Linearisation around this point gives a diagonal matrix with eigenvalues —1 in the direction („) 
and i in the direction (°). The point (0,0) is a (degenerate) source (thus unstable). Linearisation 

around it gives a diagonal matrix with the double eigenvalue 1. The point ' ^ ^ 



±) ^ 1 IS an 



unstable spiral. Indeed, linearisation gives the eigenvalues — j^ ± ^v Hi- The phase-portrait for 
system (|A2.12[) is shown in Figure [2] By using standard arguments it is easy to show that the 
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Figure 2: Phase portrait associated to (jA2.12|) showing the direction field and the eigenvectors at 
non-oscillatory critical points. The thick solid lines represent the isoclines. 



unique orbit coming out of (—1/2, 0) tangential to (_-^) ends up at the critical point I — |j — (jg) ^ 

as r — > +0O and has w < for all r. The unique orbit coming out of (—1/2,0) tangential to Cj*) 
stays unbounded in the region with dv/dw > 0. The local behaviour of these orbits give in 
the original variables give two separatrices. From the stable manifold theorem, there exist a f 
sufficiently small and a positive constant C > such that 



< e«^C for aU r < r . 



We then have 



-\<e^''C for 
2' - 



r < r 



and there are positive constants Ci and C2 such that 

6Ce^'' + Ci < z(r) <6Ce^'' + C2 for r<f 
and, in particular, there exists a constant zq such that 

lim z{r) = zq and z{r) > zq for all r <f . 



Translating this behaviour into the original variables of system (jA2.6|) . we obtain the existence 
of a unique orbit v{u) satisfying the second asymptotic behaviour in (|A2.7p . and the existence 
of a unique orbit v{u) satisfying ()A2.8p the solutions on this orbit (with w < for all r) have. 



translating the behaviour around the critical point I — |, — (:^) '' j . the behaviour (|A2.9p . 

We now infer the behaviour of orbits entering i?4 as u — > +cx3. Suppose that there are orbits 

1 , 5C 2 

with the behaviour v ~ Cu^ as |u| -> 00, then dv/du ^ 2Cu and from (|A2.6p dv/dv ' 

thus C = 1/2. To prove that all orbits, in fact, have this behaviour, we set 

v?^ M ' dz ds 



Cm^+X^ 
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and the system becomes 



dU 
ds 



-U[ V 



dV_ 
ds 



U^ + {V - 2V^ 



(A2.13) 



the transformation is vahd locahy as long a.s U ^ 0. This system has eritieal points (0, i) ,(0, 0) 
The point (O, ^) is a sink. Linearisation around it gives a diagonal matrix with 



and 



((l)^4)- 



eigenvalues — | with direction (p), and —1 with direction d') 

Linearisation around it gives a diagonal matrix with eigenvalues 1 with direction (-,). and — ^r with 



The point (0, 0) is a saddle-point. 



direction („) . Finally, the point 



is an unstable spiral. The phase-portrait for system 



(|A2.13[) is shown in Figure [31 The orbits entering the eritieal point {0,h) have thus U > and 
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Figure 3: Phase portrait associated to (jA2.13|) showing the direction field and the eigenvectors at 
non-oscillatory critical points. The thick solid lines represent the isoclines. 



V > 0, correspond to the orbits associated to solutions of (JA2.6P that at a point are in R4 with 
M > and w > for increasing u. The ones entering with U < correspond to orbits associated 
to solutions of (jA2.6[) that at a point are in i?5 with m < 0, t" > for decreasing u. Observe that 
these orbits are associated to the slow eigenvalue — | (the only orbit entering through the fast 
one, —1, is contained in the segment U ~ 0, V G {0, 1/2) and thus not correspond to solutions of 
(jA2.6p V Arguing as previously we then have that there exists a sufficiently large C > such that 
for all orbits with v{u = 0) > then 



C 
< - — r as u 

\u\ 



(A2.14) 



In particular the separatrix v{u) satisfies (|A2.14[) . In order to improve the estimate, we now take 
in (JA2.6P V = ^u^ + vi, hence vi satisfies 



S = J^+^("'-) -^*^ ^(^'-)-ifet 



(A2.15) 



Let now u > such that (u)^ > 2C where C is as above, and let L = 2C|w|5 , so that, in particular 
we can write vi{u) < LC\u\5L. Then as long as u > u and u^ < (fi) 5/(5(7) then vi{u) < Lu'^ . 
This is proved by a bootstrap argument: if vi{u) < Lus holds, then 

\R{u,v,)\<U-^ + ^lA' 
5 V|u|^ 5 / |u|5 

and hence, integrating (jA2.15|) . one gets that vi{u) < Lu's holds provided that < us < 
(«)t/(5C). . 
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Lemma 18 All orbits associated to solutions of IIA2.6\) enter R4. Those that are below the sepa- 
ratrix v come from the critical point p^ as z -^ —00. All other orbits, except for v, come from the 
region R^. 

Proof. We only need to rule out the existence of periodic orbits around the critical point. We 
then construct a spiral box invariant for decreasing z. The direction field implies that such an 
orbit can only be contained in the second quadrant. Whence, we start constructing the invariant 
set in the region R3 at a point {uo,vo) on Fi and such that < uq < Ue and vq = — -jUq. We 
continue by tracing the next point {ui,vi) on r2 vertically downwards, i.e. ui = uq, vi = ^f ^• 
The next point (u2, V2) on T^ is traced from (ui, vi) horizontally to the right;u2 = \/— 3ui = ^ , 

f2 = ui = — I — , and so on, then U4 = -^/S^^jS- and V4 = —^^^^. Observe that both uq > ^4 and 

M4 > uo hold (by the condition uq < We = ( -7? ) ) • Then, since each of the regions enclosed by 

the segment joining (u;, vi) to (u^+i, Uj+i) for i = 0, 1, 2, 3 and the isoclines have the property that 
for decreasing z an orbit can only scape the region through the isoclines, there cannot be periodic 
orbits. ■ 

The proof of the theorem follows by changing to the original variables, LemmafTBland LemmalTSl 
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